Ab initio study on the effects of transition metal doping of Mg 2 NiH 4 
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Mg 2 NiIi4 is a promising hydrogen storage material with fast (de)hydrogenation kinetics. Its 
hydrogen desorption enthalpy, however, is too large for practical applications. In this paper we 
study the effects of transition metal doping by first-principles density functional theory calculations. 
We show that the hydrogen desorption enthalpy can be reduced by ~ 0.1 eV/Kfo if one in eight 
Ni atoms is replaced by Cu or Fe. Replacing Ni by Co atoms, however, increases the hydrogen 
desorption enthalpy. We study the thermodynamic stability of the dopants in the hydrogenated 
and dehydrogenated phases. Doping with Co or Cu leads to marginally stable compounds, whereas 
doping with Fe leads to an unstable compound. The optical response of Mg 2 NiH4 is also substantially 
affected by doping. The optical gap in Mg 2 NiH4 is ~ 1.7 eV. Doping with Co, Fe or Cu leads to 
impurity bands that reduce the optical gap by up to 0.5 eV. 

PACS numbers: 71.20.-b,71.15.Nc,61.72.Bb,74.62.Dh 



I. INTRODUCTION 



The large scale application of hydrogen as a fuel de- 
pends on the development of materials that can store 
hydrogen in a dense formi Magnesium based hydrides 
are interesting candidate materials for hydrogen stor- 
age, because magnesium has a low weight. The sim- 
plest hydride, MgH 2 , has a hydrogen storage capacity 
of 7.67 wt.%. It has, however, a high hydrogen desorp- 
tion enthalpy, 0.77 eV per H2, and hence an equilibrium 
plateau pressure (10~ 7 bar) at room temperature that 
is too low for practical applications.- Moreover, the hy- 
drogen desorption/absorption reactions of MgH2 /Mg suf- 
fer from kinetic barriers, which requires operating tem- 
peratures in excess of 570 K. Various magnesium al- 
loys have been studied to improve the thermodynam- 
ics and kinetics of the hydrogen desorption/absorption 
reaction^i 4 -^^ Mg 2 NiH4, which can store 3.6 wt.% of 
hydrogen, has been suggested as a candidate material, 
primarily because of its reasonably fast hydrogen desor- 
pion/absorption kinetics i^ ' 10 : 11 ' 12 : 13 ' 14 : 15 Several theoret- 
ical studies have been dedicated to Mg 2 NiH 4 .±£»iLISii2i2fi 
However, its measured hydrogen desorption enthalpy of 
0.70 eV per H 2 r^ is barely lower than that of MgH 2 . 
It leads to an equilibrium hydrogen pressure of 1 bar 
at 510 This temperature is far too high for appli- 
cations using PEM fuel cells, for instance. Apart from 
its possible role in hydrogen storage, Mg 2 Ni is also in- 
teresting because it can act switchable mirror. 21 : 22 ' 23 : 24 
Mg 2 Ni is a metal, whereas Mg 2 NiH4 is a semiconductor 
with a band gap of 1.7-2.0 eV^^ZL This leads to a re- 
markable change in the optical properties of the material 
upon hydrogenation and dehydrogenation. Especially if 
Mg 2 Ni is applied in thin films, the optical switching can 



be fast, reversible and robust The high optical contrast 
opens up possibilities for Mg 2 Ni as a hydrogen sensing 
material 22£££L 



For both applications the reaction enthalpy of the hy- 
drogen desorption/absorption at typical operating condi- 
tions is too high. Preferable would be an equilibrium hy- 
drogen pressure of 1 bar at room temperature. To reach 
this condition a hydrogen desorption enthalpy of 0.40 eV 
per H 2 is required. A substantial amount of experimental 
work has been dedicated to study the effects of doping 
of Mg 2 NiH4 in order to reduce its hydrogen desorption 
enthalpyffi* 3 ^ 3 ^ 3 !^^^ 1 . In this context "dop- 
ing" means substituting a fairly large amount of Ni (or 
Mg) by other metals. 

In this paper we report a study on the effects of dop- 
ing Mg 2 NiH4 with transition metals by first principles 
density functional theory (DFT) calculations. We re- 
strict ourselves to the LT phase of Mg 2 NiH4, since the 
HT phase is stable only at elevated temperatures (i.e. 
T > 500K) and therefore less relevant for applications. 14 
We consider substitution of nickel by cobalt, iron, or cop- 
per in a concentration of 12.5%, which means substitut- 
ing one in eight nickel atoms. This concentration is close 
to that used in recent experiments^ 2 - Our first aim is to 
monitor the change in hydrogen desorption enthalpy and, 
in particular, to establish which dopants (if any) lead to a 
reduction of the enthalpy. We show that zero point ener- 
gies (ZPEs) associated with the hydrogen phonon modes 
considerably influence the enthalpies. Our second objec- 
tive is to study the change in optical properties that re- 
sults from doping. In particular we show that dopants in 
this concentration markedly alter the dielectric function. 
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II. COMPUTATIONAL METHODS 

First principles DFT calculations are carried out us- 
ing a plane wave basis set and the projector augmented 
wave (PAW) method ) 43 ' 44 as incorporated in the Vienna 
Ab initio Simulation Package (VASP)^^^. We use the 
PW91 generalized gradient approximation (GGA) for the 
exchange correlation functional^ The cell parameters 
are kept at the experimental values and the atomic po- 
sitions are relaxed using a conjugate gradient algorithm. 
Non-linear core corrections are applied 4^ 

To calculate accurate reaction enthalpies for reactions 
involving light elements such as hydrogen, it is essential 
to take into account the ZPE contribution. To calculate 
ZPEs we need the phonon frequencies of the materials 
involved. Phonon frequencies are calculated using a di- 
rect method^ i.e. the dynamical matrix is constructed 
from the force constants that are obtained from finite dif- 
ferences. Two opposite displacements of 0.05 A are used 
for each atomic degree of freedom. In general one needs 
to carry out such calculations on a super cell containing 
several primitive unit cells, as the force constants do not 
go to zero within a distance corresponding to a single 
unit cell. However, the unit cells of the materials stud- 
ied in this paper turn out to be sufficiently large, except 
for bulk magnesium, for which a 2 x 2 x 2 super cell is 
used. An advantage of large unit cells is that the phonon 
dispersion is small. It is therefore sufficient to calculate 
ZPEs from the phonon frequencies obtained at T. 

The dielectric functions are calculated in the or in- 
dependent particle random phase approximation taking 
into account direct transitions from occupied to unoc- 
cupied Kohn-Sham orbitals only. We neglect excitonic, 
local field and quasi-particle effects. The imaginary part 
of the macroscopic dielectric function then has the form 

^ - ^fe^S (1) 

x |(w C: k+qK,k)| 2 <5(e c ,k+q - e v ,k - hu), 

where q gives the direction of q; u,k and c, k label sin- 
gle particle states that are occupied, unoccupied in the 
ground state, respectively; e, u are the single particle en- 
ergies and the translationally invariant parts of the wave 
functions, respectively; V is the volume of the unit cell. 
Further details can be found in Ref. M. 

Almost all experimental optical data on hydrides are 
obtained from micro- or nano-crystalline samples whose 
crystallites have a significant spread in orientation. The 
most relevant quantity then is the directionally averaged 
dielectric function, i.e., e^ 2 '(ui) averaged over q. In this 
paper we only report directionally averaged dielectric 
functions. 

The Brillouin zone integrations are performed using a 
modified tetrahedron method^ All calculations on the 
hydrides use a 7 x 7 x 7 Monkhorst-Pack k-point mesh for 
sampling the Brillouin zone, and the calculations on the 
metals use a 7x7x3 Monkhorst-Pack k-point meshiS 3 - We 



use 480 bands to calculate the dielectric function. This 
number of bands includes all transitions up to 30 eV. For 
the materials containing copper a plane wave kinetic en- 
ergy cutoff of 341 eV is used, and for the other materials 
a cutoff of 337 eV. To obtain accurate formation and re- 
action enthalpies, the total energies of all final structures 
are calculated using a plane wave kinetic energy cutoff of 
700 eV. 



III. STRUCTURE AND STABILITY OF 
UNDOPED MATERIALS: Mg 2 NiH 4 , Mg 2 CoH 5 , 
Mg 2 FeH 6 , Mg 2 Ni AND ELEMENTAL METALS 

In order to assess the stability of doped Mg 2 NiH4we 
first need the total energies of the undoped hydrides and 
of all elemental metals involved. The optimized struc- 
tures of Mg 2 NiH 4 , Mg 2 CoH 5 , Mg 2 FeH 6 are given in Ta- 
ble [T] They are in good agreement with the experimental 
structures. 54 ' 55 ^ The metal atoms in LT Mg 2 NiH 4 form 
a distorted CaF 2 -type structure. Four hydrogen atoms 
are arranged around each nickel atom in a tetrahedron. 
In Mg 2 FeH 6 and Mg 2 CoH 5 the Mg and Fe/Co atoms 
form an undistorted CaF 2 -type structure. In Mg 2 FeH6 
the hydrogen atoms form regular octaeders around the 
iron atoms. In Mg 2 CoH5 the hydrogen atoms occupy five 
corners of a slightly distorted octahedron around each 
cobalt atom. 

The experimental structure of Mg 2 Ni can be found 
in Ref. The unit cell contains 12 Mg and 6 Ni 

atoms, which are basically hexagonally closed packed. 
The optimized structure given in Table U is in good 
agreement with experiment. For MgH 2 we use a previ- 
ously calculated structured For the elemental metals 
and MgCu 2 we use the experimental lattice parameters, 
i.e., a(c) = 3.21(5.21) A, a = 2.87 A, a(c) = 2.51(4.07) A, 
a = 3.52 A a = 3.61 A and a = 7.03 A for Mg, Fe, Co, Ni, 
Cu and MgCu 2 , respectively. 5 - The magnetic elements 
iron, cobalt and nickel are treated by spin-polarized cal- 
culations. The calculated pressures are small, indicating 
that it is unnecessary to explicitly optimize the lattice pa- 
rameters. We explicitly tested the latter for iron, since 
there the external pressure was largest, and obtained an 
energy gain of less than 0.01 eV. 

In order to obtain accurate enthalpies for reactions in- 
volving materials that contain hydrogen, one has to take 
ZPEs into account. All calculated total energies and 
ZPEs are given in Table HO We neglect the ZPEs of 
the elemental metals. The ZPE for magnesium is only 
0.001 eV/atom. The ZPEs iron, cobalt, nickel and cop- 
per will be even smaller, since the atomic weight of those 
elements is more than twice that of Mg. We did not 
calculate the ZPEs for Mg 2 FeH6and Mg 2 CoH5, since we 
use these compounds only to study the stability of doped 
Mg 2 NiH4 with respect to phase segregation, where the 
ZPE corrections are rather small. 

To calculate hydrogen desorption enthalpies we also 
need the total energy of the hydrogen molecule. It is cal- 
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TABLE I: Calculated atomic positions of Mg 2 NiH 4 (exp. 
Ref. [El), Mg 2 CoH 5 (exp. Ref. M), Mg 2 FeH 6 (exp. Ref. M) 
and Mg 2 Ni (exp. Ref. [13 ). In the calculations the lattice 
parameters were kept at the experimental values. 

Space group 



TABLE II: Calculated total energies and ZPEs of the undoped 
hydrides and metals per formula unit. 



Compound 


unit cell 






X 


y 


z 


Mg 2 NiH 4 


C2/c (15) 


Mg 


8/ 


0.2646 


0.4863 


0.0833 




P = 113.52° 


Mg 


4e 





0.0252 


0.2500 




a = 14.343 A 


Mg 


4e 





0.5264 


0.2500 




b = 6.4038 A 


Ni 


8/ 


0.1199 


0.2294 


0.0801 




c = 6.4830 A 


H 


8./' 


0.2088 


0.3048 


0.3041 






H 


8f 


0.1390 


0.3192 


0.8760 






H 


8/ 


0.0096 


0.2908 


0.0527 






H 


8/ 


0.1243 


0.9866 


0.0727 


Mg,CoH 5 


P4/nmm (129) 


Mg 


2a 


3/4 


1/4 







a = 4.463 A 


Mg 


2b 


3/4 


1/4 


1/2 




c = 6.593 A 


Co 


2c 


1/4 


1/4 


0.2567 






H 


2c 


1/4 


1/4 


0.4947 






H 


Sj 


0.4914 


0.4914 


0.2268 


Mg 2 FeH 6 


Fm3m (225) 


Mg 


8c 


1/4 


1/4 


1/4 




a = 6.437 A 


Fe 


4a 















H 


24e 


0.2425 








Mg 2 Ni 


P6 2 22 (180) 


Mg 


6i 


0.1639 


0.3278 







7 = 120° 


Mg 


6/ 


1/2 





0.1165 




a = 5.205 A 


Ni 


36 








1/2 




c = 13.236 A 


Ni 


3d 


1/2 





1/2 



culated using a cubic cell with sides of 13 A. We find an 
equilibrium distance of 0.7486 A, a vibrational frequency 
of 4350 cm -1 and a dissociation energy of 4.57 eV, which 
compare reasonably well with the experimental values of 
0.7461 A, 4401 cm" 1 and 4.48 eV, respectively^^ The 
0.1 eV deviation in the dissociation energy of H2 is rel- 
atively large in view of the accuracy required for calcu- 
lating hydrogen desorption enthalpies. This 0.1 eV may 
be considered as a correction to the reaction enthalpies 
discussed below. 

We calculate the ZPE for the hydrogen molecule from 
the energy levels of a Morse potential, 



E(n) 





1 






H) 




-i)] 




~ 4Dl 







(2) 



where u> is the vibration frequency and D e is the dissoci- 
ation energy. The result is given in Table HU 



IV. DOPED Mg 2 NiH 4 AND Mg 2 Ni 
A. Structure 



The unit cell of the LT phase of Mg 2 NiH 4 contains 
eight formula units. To simulate doping we replace one 
of the Ni atoms by a Fe, Co, or Cu atom, thus achieving a 
7:1 ratio between Ni and dopant atoms. In simple terms 



E (eV) ZPE (eV) 



E (eV) ZPE (eV) 



H 2 

Mg 

Fe 

Co 

Ni 

Cu 



-6.803 
-1.524 
-8.150 
-6.841 
-5.459 
-3.725 



0.266 
0.001 



MgH 2 

Mg 2 Ni 

Mg 2 FeH 6 

Mg 2 CoH 5 

Mg 2 NiH 4 

MgCu 2 



-8.983 
-9.133 
-34.511 
-29.355 
-24.053 
-9.45 



0.102 



0.852 



one can think of undoped Mg 2 NiH 4 as being constructed 
from Mg 2+ and (NiH 4 ) 4 ~ ions. The latter involve 18 
valence electrons and are closed shell ions. Upon doping 
it is likely that in the fully hydrogenated phase the closed 
shell character is maintained. This means that (NiH 4 ) 4 ~ 
is replaced by (FeH 6 ) 4 ~, (C0H5) 4 -, or (CuH 3 ) 4 -. Thus 
for a Fe atom we add two extra hydrogen atoms, for a 
Co atom one and for a Cu atom we remove one hydrogen 
atom. For all doped systems we fix the unit cell to that of 
undoped Mg 2 NiH 4 and we optimize the atomic positions. 
The external pressures on the doped systems are small, 
which indicates that the gain in energy when relaxing the 
cell volumes will not be significant. 

The geometry of the hydrogens around the Fe and Co 
dopant atoms resembles the geometry in Mg 2 FeHg and 
Mg 2 CoH5 respectively. In Mg 2 FeHg each Fe atom is in 
the center of a perfect octahedron of hydrogen atoms with 
a Fe-H distance of 1.56 A. In Fe doped Mg 2 NiH 4 the oc- 
tahedron is distorted. The H-Fe-H angles range from 80 
to 100° and the Fe-H distances range from 1.55 to 1.58 A 
in case a H atom is only bonded to a Fe atom. However, 
four of the hydrogen atoms surrounding a Fe atom also 
bond to Ni atoms, in which case the Fe-H distance is en- 
larged to 1.64-1.76 A. The hydrogen tetrahedra around 
such Ni atoms are distorted with Ni-H distances from 
1.51 to 1.80 A, whereas in undoped Mg 2 NiH 4 they are 
between 1.56 and 1.58 A. 

In the case of Co doping the distortions are much 
smaller. In Mg 2 CoHs the hydrogen atoms surrounding 
each Co atom form a four-sided pyramid with the Co 
atom just above the basal plane of the pyramid. To de- 
scribe the geometry we denoting the basal plane hydro- 
gens by Hf, and the top hydrogen by H t . The Hf,-Co-Hb 
angle is 89°, the Co-H;, distance is 1.52 A, the H&-Co- 
H t angle is 97.6° and the Co-H 4 distance is 1.59 A. The 
Co-Mg distances range from 2.75 to 2.80 A. The hydro- 
gens surrounding the Co atom in doped Mg 2 NiH 4 form a 
slightly distorted pyramid, with H^-Co-H;, angles rang- 
ing from 83.1° to 94.1° and Co-Hf, distances ranging from 
1.53 to 1.56A. The H 6 -Co-H t angle is 93.4° to 105.1° and 
the Co-H t distance is 1.57A. The Co-Mg distances vary 
from 2.69 to 2.80 A. The Ni-H bond lengths are not af- 
fected by Co doping. 

We cannot compare the geometry of the hydrogens in 
Cu doped Mg 2 NiH 4 to Mg 2 CuH3, since the latter com- 
pound is not stable with respect to decomposition into 
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TABLE III: Calculated total energies, ZPEs and hydrogen 
content of the doped hydrides, Mg2Ni 7/ / 8 TM 1 / 8 H a! , and the 
metals, Mg2Ni 7 / g TM 1 / 8 . The values are per formula unit. 



TM 

_L 1V1 


H/ ^e v f 




x \rF n ) 


Fe 


-24.952 


0.914 


4.250 


Co 


-24.677 


0.876 


4.125 


Ni 


-24.053 


0.852 


4.000 


Cu 


-23.207 


0.826 


3.875 


Fe 


-9.278 


0.100 




Co 


-9.233 


0.100 




Ni 


-9.133 


0.102 




Cu 


-8.890 


0.102 





MgH 2 and MgCu 2 ^ The Ni-H distances in Cu doped 
Mg 2 NiH4 are similar to those in undoped Mg 2 NiH4. The 
hydrogen atoms surrounding the Cu atom are located at 
three corners of a tetrahedron with the Cu atom in the 
center. The Cu-H distances, 1.62 to 1.64 A, are slightly 
larger then the Ni-H distances, 1.56 to 1.59 A. 



B. Reaction enthalpies 

In order to calculate the hydrogen desorption enthalpy 
of doped Mg 2 NiH 4 we also need the total energy of doped 
Mg 2 Ni. The unit cell of Mg 2 Ni contains six formula units 
per cell. If we replace one of the Ni atoms in this cell by a 
dopant atom, this gives a 5:1 Nhdopant ratio, instead of 
the required 7:1 ratio. We approximate the total energy 
of the 7:1 ratio by the average energy of three 5:1 doped 
unit cells and one undoped cell. All calculated total en- 
ergies and ZPEs of the doped hydrides and metals are 
given in Table [TTT1 

From the data in Table HTT1 we calculate the desorption 
enthalpy per H 2 molecule 



E, 



des 



£(H 2 ) 



2 

x 



£(Mg 2 NijTMi) 



-i?(Mg 2 NiiTMi Ha;) 



(3) 



where E(M) is the total energy of compound M and x 
is the number of hydrogen atoms in the hydride. The 
latter depends upon the dopant atom, as discussed in the 
previous section. The values of x are given in Table iLUl 
The calculated desorption enthalpy of undoped 
Mg 2 NiH 4 is 0.66 eV/H 2 without ZPE and 0.55 eV/H 2 
with ZPE. The corresponding experimental value is 
0.70 eV/H 2 jii. The agreement is fair, if we correct for the 
overestimation of the H 2 dissociation energy, mentioned 
in Sec. Mil Since the H 2 dissociation energy is overesti- 
mated by 0.1 eV, it is reasonable to assume that the des- 
orption enthalpy is overestimated by the same amount. 
This gives an desorption enthalpy with ZPE of 0.65 eV 
per H 2 , which is close to the experimental value. The 
correction is a constant shift and in the following we give 
the uncorrected results only. 



^ 0.75 r 

t5 0.7 
> 

« 0.65 



D. 



0.6 

d 

■B 0.55 
a 

% 0.5 

■B 0.45 

5 0.4 [ 
-o 



+ total energies 

x total energies with ZPE 



Fe Co Ni 

doping element 



Cu 



FIG. 1: (Color online) Desorption enthalpies -Edes (eV/Eb) of 
the doped hydrides. 



The results for Ea cs of doped Mg 2 NiH4 are given in 
Fig. m These results clearly demonstrate that the des- 
orption enthalpy can be tuned by an appropriate doping. 
The desorption enthalpy decreases considerably both for 
Fe and for Cu doping, i.e. by 84 and 71 meV per H 2 , 
respectively. However, Co doping increases the desorp- 
tion enthalpy by 28 meV per H 2 . The trends for Cu and 
Fe doping are in qualitative agreement with experimental 
data^s^ 

Fig. Q] also shows the calculated desorption enthalpies 
corrected with ZPEs. The ZPEs of all the metals are 
almost identical, see Table Hill and the ZPEs of the hy- 
drides scale linearly with amount of hydrogen atoms. 
This means that the ZPE per hydrogen atom is almost 
constant and independent of the dopant atom. There- 
fore, the ZPE correction to the desorption enthalpy per 
H 2 is 0.1 eV for all compounds studied. 



C. Stability 

Doped Mg 2 NiH4 is stable in thin films 42 In order to 
assess whether kinetics plays an important role in sta- 
bilizing these compounds, we study the thermodynamic 
stability of the doped materials with respect to phase seg- 
regation. For the dehydrogenated doped Mg 2 Ni metal we 
consider decomposition into Mg 2 Ni, bulk magnesium and 
bulk doping metal. 
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1 



Mg 2 NijTMi -Mg 2 Ni + -TM + -Mg, (4) 

8 8 o o o 

where TM = Fe, Co, or Cu. Fully hydrogenated undoped 
Mg 2 NiH4is compared to bulk nickel and MgH 2 . For the 
hydrogenated doped Mg 2 NiH4 we consider decomposi- 
tion into phase segregated Mg 2 NiH4 and Mg 2 FeHg or 
Mg 2 CoH 5 . 

Mg 2 NiiTMi He -» ^Mg 2 NiH 4 + ^Mg 2 TMH !; , (5) 

8 8 o o 



with x as in Table IIII1 and y = 5, 6 for Co, Fe, respec- 
tively. Since Mg 2 CuH 3 is unstable with respect to de- 
composition into MgH 2 and MgCu 2 , we consider for the 
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FIG. 2: (Color online) Stability with respect to phase segre- 
gation (eV / f.u.). Values are negative if the compound is 
stable, see Eqs. (gj)-© 
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hydrogenated Cu doped Mg 2 NiH4 the possible decompo- 
sition reaction 



Mg 2 Ni7CuiH a 



^Mg 2 NiH 4 + ^MgCu 2 



16 



MgH 2 



(6) 

The results are shown in Fig. [5] 

Fe doped Mg 2 Ni is thermodynamically unstable with 
respect to phase segregation into Mg 2 Ni bulk Mg and 
bulk Fe. Co doped Mg 2 Ni is a marginally unstable ma- 
terial in which segregation is favored by only ^0.01 eV 
per formula unit. Fe, Co doped Mg 2 NiH4 are thermo- 
dynamically unstable with respect to segregation into 
Mg 2 NiH 4 and Mg 2 FeH 6 , Mg 2 CoH 5 , respectively. Dop- 
ing of Mg 2 Ni with Cu leads to a stable material. Exper- 
imental work proved the stability of Mg 2 Nii^ :c Cu :c solid 
solutions) 32 ' 40 for < x < 0.85 these compounds are iso- 
structural with Mg 2 Ni. Experiment indicates that the 
hydrogenated phase decomposes into MgH 2 , MgCu 2 and 
Mg 2 NiH4i22 This is confirmed by our calculations, see 
Fig.H 

In conclusion, many of the doped phases are thermo- 
dynamically unstable. This does not need to hamper 
their usefulness, however, since kinetics plays an impor- 
tant role in stabilizing the doped compounds. The hy- 
drogen desorption temperature lies far below the temper- 
atures that is used to anneal these materials^ 2 - The Cu 
doped Mg 2 Ni metal is thermodynamically stable, and hy- 
drogenating this material can lead to a useful metastable 
compound. 



V. OPTICAL PROPERTIES 

The imaginary part of the dielectric function and the 
electronic density of states (DOS) of undoped Mg 2 NiH 4 
are shown in Figs. [3] and [4] respectively. They are in good 
agreement with the results of previous calculations ] 19 ' 20 . 
We find an indirect band gap of 1.6 eV and an optical gap 
of 1.7 eV. This is in good agreement with the experimen- 



FIG. 3: (Color online) Imaginary part of the frequency de- 
pendent dielectric functions of doped Mg 2 NiH4. 



tal optical gap of 1.7 to 2.0 e\fi&2££L The agreement is 
in fact remarkable since DFT usually underestimates the 
band gap by 30-50%. The dielectric function of Mg 2 NiIi4 
has two peaks, which can be directly related to the two 
peaks in the DOS of the valence bands. 

The dielectric functions of doped Mg 2 NiH4 are also 
shown in Fig. [3] Doping alters the dielectric function 
and, remarkably, the size of the change correlates with 
the change in hydrogen desorption enthalpy caused by 
the dopants, see Fig. [TJ This can be explained by notic- 
ing that both changes have a common cause. The changes 
in desorption enthalpy are due to changes in metal- 
hydrogen bond lengths and bond energies. Changes in 
bond energies shift the energy levels and hence can be de- 
tected in the optical spectrum. Such changes are largest 
in Fe doped Mg 2 NiH 4 . An Fe dopant atom includes two 
extra hydrogen atoms. It gives the largest perturbation 
in the Mg 2 NiH4 lattice, with almost half of the metal- 
hydrogen bond lengths being changed with respect to the 
undoped case. 

The DOS of doped Mg 2 NiH 4 is given in Figure gj To 
facilitate an internal comparison the DOS of all com- 
pounds is aligned at the bottom of the valence band. 
Besides the fundamental gap between the valence and 
conduction bands, we can identify a clear gap in the va- 
lence bands, between 2.4 and 3.6 eV below the Fermi level 
in undoped Mg 2 NiIi4. The states above this valence gap 
have a strong metal d character, whereas the lower va- 
lence states have a dominant hydrogen character. Cu and 
Fe doping introduces states in the valence gap, whereas 
all dopants introduce states in the fundamental gap. In 
the case of Cu doping these appear near the top of the 
valence band, whereas for Co and Fe doping gap states 
appear near the bottom of the conduction band. Since 
we have adjusted the amount of hydrogen upon doping, 
all doped materials are semiconducting. 

The DOS can be used to interpret the dielectric func- 
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FIG. 4: Electronic densities of state of the doped and hy- 
drides. The Fermi level is at the top of the valence bands. 
The four DOSs are aligned at the bottom of the valence band. 



tions. The decrease in the fundamental gap in the DOS 
upon doping gives a decrease in the optical gaps. The 
largest changes in the dielectric function are observed 
upon Fe doping. The highest peak decreases as compared 
to the undoped case and the valley between the two peaks 
is less deep. In addition, a distinct shoulder appears at 
low energy. Fe doping gives a clear peak in the DOS at 
the bottom of the conduction band, which yields the dis- 
tinct shoulder in the dielectric function. The two main 
peaks in the dielectric functions are not shifted upon dop- 



ing. This indicates that the dopants mainly give rise to 
additional features via the introduction of gap states, as 
can be observed in the DOS. Similar conclusions hold for 
the Cu and Co doped cases, but the perturbation of the 
Mg 2 Nifi4 DOS caused by doping is smaller than for the 
Fe doped case. 



VI. CONCLUSIONS 

Mg 2 NiH4 is a promising hydrogen storage material 
with fast (de)hydrogenation kinetics. Its hydrogen des- 
orption enthalpy, however, is too large for practical appli- 
cations. In this paper we study the effects of transition 
metal doping by first-principles density functional the- 
ory calculations. We show that the hydrogen desorption 
enthalpy can be reduced by 0.1 eV/H2 if one in eight Ni 
atoms is replaced by Cu or Fe. Replacing Ni by Co atoms, 
however, increases the hydrogen desorption enthalpy We 
study the thermodynamic stability of the dopants in the 
hydrogenated and dehydrogenated phases. Doping with 
Co or Cu leads to marginally stable compounds, whereas 
doping with Fe leads to an unstable compound. The op- 
tical response of Mg 2 NiH4 is also substantially affected 
by doping. The optical gap in Mg 2 NiH4 is ~ 1.7 eV. 
Doping with Co, Fe or Cu leads to impurity bands that 
reduce the optical gap by up to 0.5 eV. 

We study the effects of transition metal doping on the 
hydrogen desorption enthalpy and the optical properties 
of Mg 2 NiH4 by first-principles DFT calculations. The 
desorption enthalpy is reduced by 84 meV per H 2 , if one 
in eight Ni atoms is replaced by an Fe atom. Replacing 
one in eight Ni atoms by a Cu atom reduces the des- 
orption enthalpy by 71 meV/H 2 , but replacement by a 
Co atom increases it by 28 meV/H 2 . Including energy 
corrections due to the zero point motions of the atoms 
changes the absolute values of the desorption energies by 
0.1 eV/H 2 . Since however the zero point energies per 
hydrogen atom are almost independent of the compound 
studied, the relative values of the desorption energies are 
not affected. 

The thermodynamic stabilities of the doped dehydro- 
genated Mg 2 Ni and the fully hydrogenated Mg 2 NiH4 
compounds are studied by considering possible decompo- 
sition reactions. The results show that Cu doped Mg 2 Ni 
is stable, Co doped Mg 2 Ni is marginally stable, and Fe 
doped Mg 2 Ni is unstable with respect to phase separa- 
tion into Mg 2 Ni, bulk Mg and the bulk transition metal 
dopant. The doped hydrogenated Mg 2 NiH4 compounds 
are either marginally unstable in the case of Co or Cu 
doping, or, in the case of Fe doping, clearly unstable. 
Kinetic barriers could be be sufficiently high to stabilize 
metastable doped compounds since the hydrogen desorp- 
tion temperatures are smaller than the temperatures used 
to anneal these materials. Nevertheless, thermodynam- 
ics indicates that Cu is the most promising candidate to 
lower the hydrogen desorption enthalpy of Mg 2 NiH4. 

By calculating the dielectric function within the ran- 
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dom phase approximation we study the effects of doping 
on the optical properties of Mg 2 NiH 4 . The changes in 
the dielectric function can be interpreted in terms of the 
electronic densities of states of the corresponding com- 
pounds. The dopant atoms introduce states in the fun- 
damental gap, as well as below the valence Ni d-band. 
These states cause a shift in the onset of absorption to 
lower energy by up to 0.5 cV and they decrease the rel- 
ative heights of the peaks in the Mg 2 NiH4 absorption 
spectrum. The sizes of these changes correlate with the 
change in the hydrogen desorption enthalpy caused by 
the dopants. Fc doping causes the largest disruption in 
the Mg 2 NiH 4 lattice, and the largest change in the opti- 
cal properties. 
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